Collective excitations and the gap in spectrum of the degenerated 

Bose gas 

V.B. Bobrov, S.A. Trigger 
Joint Institute for High Temperatures, 
Russian Academy of Sciences, 
13/19, Izhorskaia Str., 
Moscow 125412, Russia; 
email: satron@mail.ru 

Model of the degenerated weakly non-ideal Bose gas is considered without sugges- 
tion on C-number representation of the creation and annihilation operators with zero 
momentum. The "density-density" correlation function and the one-particle Green 
function are calculated on basis of the suggestion about C-number representation of 
the operator of particle density in the Bose condensate. It is shown, that the pole 
in the "density-density" Green function determines the Bogolyubov's spectrum of 
the collective excitations, which have the phonon-roton form. At the same time the 
spectrum of the one-particle excitations possesses a gap, whose value is connected 
with the density of particles in the Bose condensate. 



PACS number(s): 05.30.Jp, 03.75.Kk, 03.75.Nt, 05.70.Fh 



I. INTRODUCTION 



The fundamental ideas and approach on explanation of superfluidity in Hell have been 
formulated by Landau [1,2]. The crucial role in these papers play the statement about the 
phonon-roton spectrum of the elementary excitations and the Landau condition for super- 
fluidity. The microscopic model of the weakly non-ideal Bose gas below the condensation 
temperature T has been constructed by Bogolyubov [3] by use the concept of the C-number 
character of the creation Oq and annihilation ao operators and diaganalization of the ap- 
proximate Hamiltonian, executed on basis of this concept. 
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The spectrum of the elementary excitations Hcu(q) in this model corresponds for a small 
q to the real phono n-roton spectrum in Hell at T < T , which follows from the neutron 
scattering experiments [4,5]. For large values of q this spectrum corresponds to one for 
the free particles e q = h 2 q 2 /2m. The respective diagram technic for the Bose systems at 
T < To has been developed by Belyaev [6]. It is necessary to stress that in the early papers 
Landau and also Bogolyubov mentioned possibility of the existence of the spectrum with 
a gap, but later they omitted this idea, because contradiction one with the observations of 
the phonon branch of spectrum at small q in the neutron scattering experiments. Recently 
the arguments and the preliminary estimations on coexistence of two branches of spectrum 
- the Bogolyubov's collective, and a new - one-particle with a gap, have been done in [7-10]. 
The proof of this concept can be done by the various approach, including the dielectric 
formalism, developed for the Bose system at T < T in [11-13]. In the present paper this 
statement is justified in the framework of the hypothesis about C-number behavior of the 
density of the particles in the state with p = 0. 



II. DENSITY-DENSITY CORRELATION FUNCTION AND MODIFIED 

C-NUMBER APPROACH 

Let us consider the weakly non-ideal Bose gas at T < T . Then, the distribution function 
in momentum space 

f(p) = H a p) (!) 

by analogy with the case of an ideal Bose gas takes the form 

/(p)HW,o + /J(i-W, (2) 

Here a+ (a p ) are the operators of creation (annihilation) of a particle with the momentum hp, 
the angle brackets mean averaging on the grand canonical ensemble with the Hamiltonian 
H 

H = J2 £ P a P a P + ^ U ^ fl p + i- 9 /2 a J 2 +«,/2 fl P2-9/2flpi+9/2 (3) 

p q,pi,P2 

Here \i is the chemical potential, the function u(q) is the Fourier-component of the interaction 
potential between the particles, N = J2 P a p a pi N = a^ao - the operator of the number of 
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particles with the momentum equal zero - (the "condensate"), /J = (a+ a p ) - the one-particle 
distribution function for the particles with the non-zero momentum (the "over-condensed" 
particles). 

Let us consider the temperature "density- density" Green function x(g, if2) (q ^ 0) 

1 

X(q,ity = y((p q \p- q ))«i (4) 

where p q = S p Op_ q / 2 a p+q/2 is the Fourier-component of the particle density operator. Then 
the function x(g, iQ) can be represented in the form 

X(q, ifi) = -y- S p F(p, q, iQ), F(p, q, iO) = (( a+ q/2 a p+q/2 | p_ q )) in (5) 

The equation of motion for the function F(jp, q, in) with the Hamiltonian H (Eq. (J3j)) reads 



(in + £p- q / 2 - ^p+q/2) ^(P, q, ity = /p-q/2 ~ /p+q/2 ~ 

y ( a i+k-q/2 a J 1 -k/2 a Pl+k/2a P +q/2 ~ ^-q^^-k^Pl+k^p-k+q^) | P-q))id§) 

k pi 

Taking into account Eq. (j2J), from Eq. (J5j) we find that the function F(p, q, ifl) has the sin- 
gularities at p = ±q/2. Therefore, the density-density function x(g, in) can be represented 
in the form 

X (q, in) = 1 F(q/2, q, zfi) + i F(-q/2, q, ifi) + 1 ^ F T (p, q, ifi) (7) 

p^±g/2 

The index T means, that the respective function concerns to the "over-condensed" particles. 
Then in the last term in Eq. (J7]) we can change summation by integration on momenta. The 
functions F(±q/2, q, in) extracted above satisfy, according to Eqs. (12), (El), to the equations 
of motion 

(in - e q ) F(q/2, q, in) = (N ) - — ^ u(k) (a^a+ _ k/2 a 6 / pi+ k/2a q - 

fc^O pi 

a a p!-k/2 a Pi+k/2a q -k) I P-q))i(l (8) 

(ifi + £ g ) F(-q/2, q, ifi) = -(N ) U< ^> X^ ( a k-q a P " 1 -k/2 a Pi+k/2a Q - 

a_ q a pi _ k /2 a pi+k/2a-k) I P-q))iSl- (9) 
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Let us consider further the case of strongly degenerated gas, where T — > 0. In this case we 
accept (N ) — > (N). According to the Bogolyubov's procedure, let us extract in the right 
sides of Eqs. (jHJ),© the main terms, which are determined by the maximum quantity of the 
operators and ao- In the limit of a strong degeneration and a weak interaction the other 
terms can be omitted. Then for q 7^ Eqs. (JEJ),© read 

(ittTtq) F(±q/2,q,ift) = ± ((N ) + ^u(g)(( (a^a q a + a^ai q a Q a ) | P- q ))«ij (10) 

Let us suppose that the operator of the quantity of particles in the "condensate" No is 
the C-number (No =< N >). At the same time the operators ao are not the C-numbers 
(aoCLQ — a^ao = 1) in contrast with the theory of Bogolyubov. In this case from Eqs. ffTU]) 
straightly follows 

(iQ =F e q ) F(±q/2, q, in) = ± ((No) + ^u(q){F(q/2, q, in) + F(-q/2, q, iQ}^ (11) 
From Eqs. (fTTj) one find the solutions for the functions F(q/2,q,in) and F(—q/2,q,iQ) 

w«> ~ m^r Fi ~ M - - tw-Swr (12) 



hu(q) = ^Je 2 q + 2n u(q)e q (13) 

where uq = (No) /V is the average density of particles in the "condensate". The relation 
( JTBl for the spectrum exactly corresponds to the known Bogolyubov's expression. Inserting 
( TT2j) in ffTTl) and taking into account that for a strong degeneration and a weak interaction 
the contribution of the functions F T (p, q, in) is negligible, we arrive at the expression for 
the "density- density" Green function [11,12] 

* ift ) = PFMf < 14 > 

The function x(g, in) can be continued analytically to the upper semi-plane of the complex 
z, where this function coincides with the retarded "density-density" Green function x R (Qi z ) 
(with the change in on Tiz). As it is well known, the singularities of the function x R (g, z) 
determine the spectrum of collective excitations in the system. Therefore, in the suggestion 
about C-number behavior of the operator N we reproduce the Bogolyubov's result for the 
spectrum of excitations in the degenerated and weakly interacted Bose gas. 
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III. ONE-PARTICLE SPECTRUM WITH THE GAP 

The necessary conditions for correctness of the above derivation are the conditions 
lirriT^o{No) = (N) and finiteness of the function f£ for p —* 0. Let us show that these 
conditions are fulfilled in the degenerated weakly non-ideal Bose gas. For this purpose we 
calculate the one-particle Green function g(p,iQ) for p ^ 

g(p,iQ) = ((ctp | a p ))iQ (15) 

The equation of motion for the Green function g(p, ifl) for p ^ reads 

(itt-Ep + fj,) g{p,itt) = 1 + TT^2u(k)^2((a+ i+k a bfpi a p+k | a+))i 



V , 

k pi 



I HI 



1 + < y > u(0)g(p,iQ) + ^p^2u(k) ^(( Qp 1+k a Pl a P +k I a+)) iQ . (16) 

fc^O Pi 

where u(Q) = u(p = 0) = lim p ^Qu(p). As for the "density-density" response function, we 
consider the case of a strong degeneration and a weak interaction and extract in the right 
side of Eq. (fl6|) the terms with maximum quantity of the operators Oq and ao- Then from 
Eq. ([16]) one find 

(it! - e p + /i (0) ) g{p,iQ) = 1 + yU(p)((a^a p a | aj)) in (l - S p>0 ), (17) 

where fx^ = /i — mt(0) and n = (N) /V is the average density of the particles. Further, 
suggesting again the C-number behavior of the operator No and taking into account that in 
the approximation under consideration jj^ = 0, we find for p ^ 

9 < p - ffl > = ^- < 18 > 

The expression for the spectrum E p is equal 

E p = e p + n u(p), (19) 



From Eqs. (fl8|h(fI9j) for T < U directly follows 

/pT= exp(^ p )-l" (20) 
where (3 = 1/T. Therefore, the function /J is finite for p — > 0. Moreover, 

f v - 0, (21) 

for the limit p — > 0, T —>■ in opposite to the case of the ideal gas, where f p ,ld — > oo for 
p — > 0. Accordingly, the relation ([2]) is valid for the interacted system. Besides that, in the 
one-particle excitation spectrum appears the gap A = E p ^ = n u(0). 
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IV. CONCLUSIONS 

Therefore, even a weak interaction between Bose particles leads to the drastic distinction 
with the case of the ideal Bose gas. This distinction takes place not only for the collec- 
tive excitations, which are described by the "density-density" correlation function, but also 
in the distribution function and the one-particle excitations of the over-condensed parti- 
cles. In the case under consideration, for strong degeneration T <C T , the influence of the 
over-condensed particles and the temperature effects are small. However, the conclusion 
about existence of the gap is evidently valid in the region of temperatures T < T . In this 
connection we have underline that the spectrum E p (fl9l) satisfies the Landau criterion of 
superfluidity for transitions between the condensed and over-condensed particles. For the 
transitions between the over-condensed particles this criterion violates, as it was mentioned 
by Bogolyubov. 
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